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We present the status of the CKM matrix parameters in the framework of the Standard Model. 
We perform a model independent analysis to set constraints on additional effective parameters 
accounting for possible New Physics effects and to evaluate the present allowed space for these 
effects both from Bd and Bs mesons. 



The unitary Cabibbo-Kobayashi-Maskawa (CKM) matrix|I^ describes the mixing of the quark 
flavors within the framework of the Standard Model (SM). Profs. Kobayashi and Maskawa 
have just been awarded the Nobel prize for their early 70's work on such a 3 x 3 (3 quark 
generations) unitary matrix that accounts for violation of CP symmetry through Electro- Weak 
(EW) couplings. It has 4 real parameters, amo ng w hich one single non- vanishing phase. We 
employ an exact Wolfenstein-like parametrizationl2EI that describes the strong hierarchy in these 
couplings where unitarity holds to an arbitrary power of the Cabibbo angle A = sm{9c), it is 
also re-phasing invariant: 



The parameter A is accurately determined (at 0.3 % level) from super-allowed nuclear transitions 
(|K(d|) and in semi-leptonic kaon decays (|T4s|). The parameter A {\Vcb\) is measured from 
charmed B semi-leptonic decays with an accuracy at the level of 3 %. The apex of the Unitary 
Triangle (UT), i.e. the complex number (p + if]), is less constrained. 

The accurate measurement of these parameters and especially of the UT sides and angles, 
possibly in a redundant way, allows to check the consistency of the Kobayashi-Maskawa (KM) 
mechanism within the SM. Any significant departure could suggest contributions from New 
Physics (NP). The challenge, both for experimentalists and theorists, is that precise extraction of 
observables related to these EW parameters is complicated by the presence of strong interaction 
effects. 

We perform a global fit to the CKM parameters within a frequentist approach including a 
specific treatment to deal with theoretical uncertainties (i.e. flat likelihood d la Rfit^, where 
we only use the observables from K and B sectors on which we have a good theoretical control, 
to avoid to claim pseudo departures from SM. Table [T] displays the various key ingredients used 
(more details on the world averages (WA) exp. and theo. inputs and related references are given 
atEl). Among all these observables, only the branching ratio (BR) of the B"^ P^P^ channel 
updated by the BaBar collaboration^ is a new input since our last summer 2008 update. 

Several hadronic inputs are mandatory for the fits. They mainly limit the precision on 
the determination of the observables involving processes with loops such as Am^, Am^, \£k\, 



Phys. params. 


Expcrim. input 


Theory method/ingredient 


\Vcb\ 
\Vub\sL 

\yub\lept. 

Arris 
Arrr^ 


super-allowed /3 decays 
/Cjs SL kaon decays (WA, FlavianetEl) 
B Xclu (HFAgEI excl. + incl.) 
B Xulu (HFAgEI excl. + incl.) 

BR(i3+ — > T^v) annihilation (B-factories) 
Bs ~ Bg mixing (CDF II) 
- mixing (HFAgEI) 


Towner and Hardy (08) 
f^^(0) = 0.964(5) (RBC-UKQCD (07)) 
40.59(38) (58) x lO'^ (FF and/or OPE) 
3.87(9)(46) X 10-3 (FF and/or OPE) 
and own syst. treatment 
[fB, = 228(3)(17) MeV, /s,//b, = 1.196(8)(23)] 
Ai?=2 amp. [Bg = 1.23(3)(5), /s^ , rfit , rjg] 
AB=2 amp. [BbJBb^ = 1.05(2)(5), 

fsJfBa, Vb] 


7/03 


KK mixing (PDG 08''': KLOE, NA48^TeV...) 
Charmonium B decays (HFAGEI) 
B tttTjPp, pTT (B-factories; rates + asym.) 
B- ^ (B-factories: rates -|- asym.) 


AS=2 amp. [Bk = 0.721(5)(40), rjcc , Vet, m] 

Isopsin SU{2) (Gronau, London (90)) 
GLW/ADS/GGSZ 



Table 1: Various relevant inputs to the CKMfitter global fit. Many LQCD inputs in these table are from our own 
average (see text). The upper (lower) part of the table corresponds to CP conserving (violating) parameters. 



and also the tree decay ^ t^v. The hadronic contributions to Ki^ decay are surprisingly 
under excellent control. We mostly rely on lattice QCD (LQCD) simulations to estimate these 
quantities, since the accuracy of such first-principle computations can be improved in a controlled 
way (at least in principle). The presence of results from different collaborations with various 
statistics and systematics makes it all the more necessary to combine them in a careful and 
reproducible way. It has been pointed ouP that "i/ experts cannot agree, it is unlikely the rest 
of the community would believe a claim of NP". Therefore we have recently set up our own 
average of these results §. 

Figured] (Left) shows the global CKM fit results in the (/),?/) plane. The CKM parameters 
are: A = OMWtomAi, ^ = 0.22521 ± 0.00082, p = 0.139t|]:[}27, and fj = OMltoml. A good 
overall consistency at 95 % CL is seen, probing the fact that the KM mechanism is at work for 
CP violation and dominant in B decays. It is also visible that there is a tension between the 
measurement of sin(2/3) from charmonium B decays and the determination of \Vub\ from the 
decay B^ — > t~^v. When removing one of the last paramaters from the global fit, the at 
minimum drops respectively by 2.3 and 2.4 o". 

" We apply the averaging procedure*^ 

• First of all, we collect the relevant calculations of the quantity that we are interested in: we take only 
unquenched results with 2 or 2-1-1 dynamical fermions, from published papers or proceedings. In these 
results, we separate the error estimates into a Gaussian part and a flat part (Rfit). The Gaussian part should 
collect the uncertainties from purely statistical origin, but also the systematics that can be controlled and 
treated in a similar way (e.g., interpolation or fitting in some cases). The remaining systematics constitute 
the Rfit error. If there are several sources of error in the Rfit category, we add them linearly, keeping 
in mind that in many papers in the literature, this combination is done in quadrature and the splitting 
between different sources is not published. If Rfit is taken stricto sensu and the individual likelihoods are 
combined in the usual way (by multiplication), the final uncertainty can be underestimated, in particular 
in the case of marginally compatible values. 

• We correct this effect by adopting the following averaging recipe. We first combine the Gaussian uncer- 
tainties by combining likelihoods restricted to their Gaussian part. Then we assign to this combination 
the smallest of the individual Rfit uncertainties. The underlying idea is twofold: (1) the present state of 
art cannot allow us to reach a better theoretical accuracy than the best of all estimates, and (2) this best 
estimate should not be penalized by less precise methods (as it would happen be the case if one would 
take the dispersion of the individual central values as a guess of the combined theoretical uncertainty). It 
should be stressed that the concept of a theoretical uncertainty is ill-defined, and the combination of them 
even more. Thus our approach is only one among the alternatives that can be found in the literature. In 
contrast to some of the latter, ours is algorithmic and can be reproduced. We found a very good agreement 
between our previous inputs (taken from lattice reviews) and our current set (obtained from the above 
recipe). 



-0.4 -0.2 0.0 0.2 0.4 0.6 0.8 1.0 20 40 60 80 100 120 140 160 180 

P a (deg) 

Figure 1: 95 % CL individual and global constraints in the {p,fj) plane from the global CKM fit (Left). The red 
hashed region of the global combination corresponds to 68 % CL. CL profile for a with the present world average 

of the 3 B ^ Tvn, pp, pn channels (Right). 



This t ensi on is mainly originated from the recent BR(i?^ t^v) measurements by BaBar 
and Belle^2IHl_ pj[ these measurements are consistent and their WA is (1.73 ib 0.35) x 10~^, 
while our global CKM fit predicts it to be at a lower value of (O.SOlQQg) x 10^^. Such a higher 
BR is not necessarily accommodated for by models with 2 Higgs boson 

doubletstini (2HDM). 

In addition one can see on Fig. [1] that both semi-leptonic and purely leptonic B decays \Vub\ 
determinations agree pretty well. In between the 2 semi-leptonic methods sin(2/?) prefers the 
exclusive one, while the inclusive one is still compatible in the CKMfitter approach. Doing the 
computation of the ratio of the BR of this B annihilation decay over the mixing parameter Amc^ 
removes the dependance to the decay factor /s^. The combination of these 2 constraints releases 
therefore partially some LQCD related uncertainties and gives a direct access to the parameter 
Bb^- When doing so we obtain the value Bb^ = 1-18 it 0.14 that is 2.7 a away from the CKM 
global fit: 0.521°; Jf. The tension arising from the BR(i3+ —> t^u) is clearly not yet an evidence 
for NP, but it motivates more accurate measurements at BaBar and Belle and at possible future 
super-B factories. 

It has been suggested^^ that the recent LQCD improvements in the determination of the 
parameter Bk alights a so far neglected additional multiplicative factor in the determination of 
the parameter \eK\-, this is the so called Ke parameter computed and estimated to be equal to 
0.92 lb 0.02 . This factor accounts for CP violation effects in K — K mixing and may hint for CP 
violation contributions originated from NP. The computed value of |ei<-| from this recent work and 
within the SM is (1.78±0.25) x 10"^ while the current experimental WA^lis (2.229±0.10) x lO'^. 
This suggests an additional tension at the level of 2 o" mainly with respect to sin(2/3). Our fiP^, 
even while accounting for k^, shows that the uncertainty of |ei^| is rather likely to be of the order 
of 0.5 X 10"'^. This tension arises while dealing with pure convoluted Gaussian uncertainties for 
all the parameters and including all the uncertainties on LQCD computations, that are obviously 
not overwhelmed by statistical effects. It therefore vanishes while using the Rfit procedure. 

Figure [T] (Right) shows that the angle a is now determined with a good accuracy, at the 
level of 5 % or less: a = (89.0l42)°, while the angle /? is measured within a precision of 4 %. 
The isospin analysis on the pp channels almost fully drives it. It is in excellent agreement with 
the global fit (95.6lgg)° (without the related measurement in the fit) and the uncertainties 
have been reduced by more than 20 % with respect to last summer. This is due to the new 
measurement on the BR(i?^ P^P^) by BaBar'^ that dominates the WA for this observable. It 
has increased from (18.2 it 3.0) x 10~^ up to (24.0 it 1.9) x 10~^. In the pp system, the Penguin 



to Tree amp litude ratio is much more favorable than in the case of charmless B decays to pvr 
and 7r7r'2112J^ allowing therefore a relatively smaller |Aa| isospin bound. 

The BR of both channels p~^p^ and p~^p~ are now very similar^l and almost 25 times as 
big as that of p^p^ (the Penguin transition), the B and B related isospin amplitudes triangles 
are basically flat and do not close, i.e. for B : \A'^^\/\/2 + < |^+°| (but this is stih 

consistent within uncertainties). As a consequence the mirror solutions that possibly arise while 
experimentally measuring the effective angle cxeff (Penguin dilution), are degenerated into a 
single peak. As it can be seen on Fig. [1] the expected 8-fold ambiguities from the isospin analysis 
degenerate into the only 4 Aa geometric solutions, in the vicinity of 0°, it 90°, and it 180°. 

The isospin analysis for the pp system is performed using^i the 3 BRs, time-dependent CP- 
asymmetry parameters C"' , , C^^, and C^^, and the 3 longitudinal fractions {/l) of these 
VV channels that are not stricto-sensu CP-eigenstates, thought the /l are very close to 1 which 
eases the analysis. The a angle is determined to be (89.9 it 5.4)° and the isospin bound Aa 
close to 0° with a good accuracy: (1.4 it 3.7)° (at summer time we had: a = {90.9~^fjg)°). To 
test what is the expected uncertainties for this measurement, we have performed 1000 pseudo 
experiments (toys). We have generated the above experimental observables with ±1 a around 
their best fitted value (from global fit), where the a are the currently measured uncertainties. 
We measure that the average expected uncertainty is 7.5°, slightly higher than the 5.4° that 
we measure. The uncertainty distribution has a long tail up to about 20°, it corresponds to 
revival of pseudo mirror solutions, above the 1 a CL(q). About 34 % of the toys where isospin 
triangles close and have similar uncertainties or higher than that of last summer configuration. 
This is a message for future experiments, such as LHCb, that better uncertainties of the various 
pp observables may not necessarily lead to better accuracy on a. 

Due to the reached precision, it is legitimate to investigate for possible isospin breaking 
effects^^ beyond the Gronau-London SU{2) method. Not all the breaking effects can be cal- 
culated at present, but we can list a few of them: the u and d quarks have different electric 
charges and masses (breaking of the order: (m„ — vrid) / ^qcd ~ 1 %), the isospin transitions 
A/ = 5/2 may be no more negligible, we may need to extend the basis of EW- Penguin operators: 
Q7,...,io (Aa£;vKP ~ 1-5°), the mass and isospin eigenstates are different {p — io mixing at the 
level of 2 %), the p natural width is large enough such that 1 = 1 contributions are possible 
(©(Fp/mp) ~ 4 %) ... There are possible ways out such as exploiting the B^ — > K*p^ channels 
through SU(3) constraints. In order to break the triangle closure we apply the procedure as 
described iiv^. The amplitudes ^4"'"'' and A~^^ are corrected by additional Tree (Ay) and Penguin 
(Ap) contributions weig hted as: V2AA+^ = VudV*f,ATT+~ + VtdV^lApP+- (the strong phases 
are set arbitrarily). We tested |A^"'"''| as big as 4, 10, and 15 %. The 2 first corrections break 
SU{2) at 90° and restore it in the vicinity of 0°, while the largest is needed to restore it at 
the a SM solution. Anyway when combining the tttt and pvr the determination on a is mostly 
unaffected at 1 cr CL. 

We have updated the constraint on |Vtd/Vts| accessible through the ratio of branching 
ratio for B — > decays, where V holds respectively for {p,uj) and K* vector mesons. These 
penguins processes complement the box diagrams involved in the measurement of Am(^d,s)- Any 
inconsistency in between the 2 approaches would teach us in which direction to look for NP. We 
use the parametrization for hadronic effects as described in^. The sophisticated description of 
the amplitudes has non trivial sensitivity to the CKM parameters. Our new analysis benefits 
from the recent updated BR measurements of all of the above decays^. The improvement is 
such that at 95 % CL these new measurements constrain the {p,fi) plane as accurately as Am^ 
alone, and at 68 % CL they have similar precision as that from Am(^d,s) 95 % CL. 

There has been a standing issue due to apparently non SM BR measurements for leptonic 
decays of Dg mesons^^^^^^, by the B-factory and the CLEO-c experiments. These decays give 
access to the measurement of the decay parameter fj^g and to \Vcs\- The charm sector, where 



nic ~ Aqcd, is an ideal laboratory to validate LQCD against experiment. The recent most 
accurate BR measurements by CLEO-cff^on annihilation decays Dg — > (t, allow to compute 
fDs = (259.5 ± 6.6 ± 3.1) MeV, while our average on LQCD results is (246.3 ± 1.2 ± 5.3) MeV. 
There is still some discrepancy at the 2 a level, but it is almost twice as less as what it used to 
be. Converting this into a \ Vcs\ determination and averaging CLEO-c and LQCD measurements 
of /ds, one computes \Vcs\ = 1.027 ± 0.051, in good agreement with the global fit that yields 
0.97347 lb 0.00019. This comparison alighted a 2 o" tension one year ago and the measurements 
led to a unitarity violation of the CKM matrix 13]. 

We also updated the constraint from the measured BR of the tt^vv rare decay, 

for which a recent update of the E789 and E949 experiments has been done with 5 signal 
candidates . We parameterize the BR using the calculations by Brod and Gorbahn at NLO 
QED-QCD and accounting for EW corrections to the charm quark contribution. The global fit 
predicts BR=(0.811 



+0.027 
-0.021 



exp. 



lb 0.096t/ieo.) X 10 while the experiments measure {l.Ti_i 



05 J 



-|^g-io_ 'pj^g agreement is good and the constraint in the (/),?/) plane is such that in the vicinity 
of the point (1,0) a non negligible area is forbidden at 95 % CL for the first time. This effect 
clearly motivates a O(IOO) signal event experiment, such as the future NA62. 

Finally we reiterat J^^^ the analysis to compute the constraints set on NP from 5^=^ .^-meson 
mixing. We consider that NP only affects the short distance part of the /S.B = 2 transitions. 
In addition we assume that the tree-level mediated decays proceeding through a Four Flavor 
Change get only SM contributions (SM4FC hypothesis: h — > qiqjqk (^ / J / ^))! the observables 
\Vij\ (including — > r+i/), 7, and 7(a) = vr — Pec — ol are not affected by the NP contribution 
and can be used in a (SM+NP) global fit to fix the SM CKM parameters. We also consider 
only 3 generations of quarks. The oscillation parameters, the weak phases, the semi-leptonic 
asymmetries and the i?-meson lifetime differences are affected by the phase and/or the amplitude 
of the NP contribution and allow to constrain the NP deviation to SM quantified through out 



the model- independent parametrization: (i?q|M^2 ^ l-^i 
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Figure 2: 68 % CL contours for A^'' in Bd — Bd system (Left) and in Bs — Bs system (Right). 



In Figure [2] we present the deviations to the SM (A^^ = 1) using the intuitive Cartesian 
coordinates parametrization A^^ = (Re + i Im) A^^. This parametrization is statistically 
more robust as uncertainties have Gaussian behavior in the vicinity of |Aq^| = 0. In the B^ 
case, the tension in between sin(2/3) and iVufelri^ pushes the best fitted A^^ 2.1 a away from 
the SM point (while it is only 0.6 a away when B~^ — > t~^v is removed). In the case of Bg, 
the deviation is 1.9 cr, it's mainly driven by the recent TeVatron measurements of (2/?^, AFg)^. 



This measurement is performed with the time dependent analysis of the decay Bg J /tptj). It 
deviates by 2.2 a from the SM expected value. In both cases Amq=d,s constrain the modulus 
I A^£^ ^1 to be in the vicinity of 1 or below, this is the evidence of the KM mechanism dominance 
for the sensitivity to NP effects. If one tests the Minimal Flavor Violation (MFV) scenario (i.e. 
no additional NP phase and Yukawa couplings only: ImA^^ = and A^^ = A^^), no tension 
with respect to SM is observed, as theses effects arise at the present time only through EW 
phases: sin(2/3) vs. |V^fe|r!/ and (j)s, in both Bq=ci,s systems. 

To conclude the KM mechanism is at work and dominates the sensitivity to CP violation 
and to NP in the b quark sector. Anyway there is still substantial room for NP both in 3^- 
meson and S^-meson physics. Some few deviations to the SM global fit exist at the present time 
and at most at the 2 a level. It is therefore fundamental to finalize the analyzes of the present 
B-factory datasets and to wait for the next generation experiments at the LHC (huge h quark 
cross-section production), or at the future super-B factories, at KEK and possibly at Frascati 
{C = 10'^^~'^^cm~-^s~^). They will allow for high precision measurements of rare effects. Finally, 
continuous progress in LQCD are currently achieved, but even more accurate calculations, in a 
coherent motion of that community, are mandatory and expected to fully exploit the potential 
of the physics program in that field. 
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